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! ABSTRACT 

Context. The warps in many spiral galaxies are now known to asymmetric. Recent sensitive observations have revealed that 
, asymmetry of warps may be the norm rather than exception. However there exists no generic mechanism to generate these 
^ asymmetries in warps. 

Aims. To provide an explanation for the generation of asymmetric warps in disk galaxies 
. Methods. We have derived the dispersion relation in a compact form for the S-shaped warps(described by the m = 1 mode) 
00 ' and the bowl-shaped distribution(described by the m = mode) in galactic disk embedded in a dark matter halo. We then 
' performed the numerical modal analysis and used the linear and time-dependent superposition principle to generate asymmetric 
warps in the disk. 

Results. On doing the modal analysis we find the frequency of the m = mode is much larger than that of the m — 1 mode. 
The linear and time-dependent superposition of these modes with their unmodulated amplitudes (that is, the coefficients of 
' superposition being unity) results in an asymmetry in warps of ~ 20 — 40%, whereas a smaller coefficient for the m = mode 
Oh results in a smaller asymmetry. The resulting values agree well with the recent observations. We study the dependence of the 
^ , asymmetry index on the dark matter halo parameters. This approach can also naturally produce U-shaped warps and L-shaped 
' warps. 

' Conclusions. We show that a rich variety of possible asymmetries in the z-distribution of the spiral galaxies can naturally arise 
d , due to a dynamical wave interference between the first two bending modes(i.e. m = and m = 1) in the disk. This is a simple 
^ ' but general method for generating asymmetric warps that is independent of how the individual modes arise in the disk. 
• I— I 
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1. Introduction generate an S-shaped warp in the galactic disk (Battaner 

et al. 1990). But a very high value of the magnetic field 
The disks in most spiral galaxies are not flat and their ^^^^.^^^ p^^^^^^ ^^^p ^^^^ 4_g ^.^^ scalelengths. 

outer parts often show warping of the material away ai..-i j.- u j. j. 

^ ^ ° Alternatively gas accretion could generate warps repeat- 

from the galactic midplane. Warps are mostly seen in the ^^j^^ exa.mp\e seen in the cosmological simulations 

atomic hydrogen gas (HI) in the outer parts of a galactic ^.^^ formation (Scmclin & Combes 2005). 

disk (Sancisi 1976, Bosma 1978), but there is also a strong 

observational evidence for optical warps in old stars, as in One plausible mechanism to get long-lived warps is 

NGC 4565 (vand der Kruit & Searle 1981). A statistical to use the modal approach first proposed and studied 

analysis of a large number of edge-on spirals has shown by Hunter & Toomre (1969). Sparke & Casertano (1988) 

that optical warps are seen in a significant fraction (about (hereafter, SC88) followed Hunter & Toomre's (1969) ap- 

40 %) of the spiral galaxy disks (Sanchez-Saavedra et al. proach but also included a dark matter halo and achieved 

1990, Reshetnikov & Combes 1998, 1999). the above goal of avoiding differential precession. They 

Thus the warps are a common phenomenon and hence ^ave shown that in the presence of an oblate halo po- 

must be either long-lived or be excited repeatedly. Tidal tential a system of concentric circular rings, representing 

interactions have been proposed as a triggering mecha- the galactic disk, is able to accomplish a configuration 

nism for warps, however since many isolated spirals- such ^here all the rings are synchronized to precess with a uni- 

as NGC 4565 also show warps (Sancisi 1976), this cannot ^rm precession frequency. Such modes are called a normal 

be the main triggering mechanism for warps. Anisotropic modes, and in this case since the halo is included, it is a 

pressure applied by the intergalactic magnetic field can modified-tilt mode of the disk because it is derived by 

modifying the trivial tilt mode solution of the linearized 

Send offprint requests to: K. Saha equations of motion of an isolated system of rings. 
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The warps were initially deduced to have an S-shaped 
distribution, with equal amplitudes on both sides of the 
mid-plane. This makes their theoretical analysis easier. It 
has long been known, however, that quite a few galaxies 
show asymmetric warps so that the amplitudes on the 
two sides are appreciably different. For example, our own 
Galaxy has an asymmetric warp: the mean plane bends 
up to 4 kpc above the plane of the inner Galaxy on the 
northern side of the disk, while it reaches down only half 
on the southern side and then turns back up toward the 
inner disk plane (sec Burton 1988). More recent, sensitive 
observations have revealed that asymmetry of warps may 
be the norm rather than the exception as shown by the 
HI gas and optical R-band observations (Garcia-Ruiz et al. 
2002). A detailed catalog of the measure of asymmetries 
in optical warps seen in edge-on galaxies can be found 
in Sanchez-Saavedra et al. (2003). Theoretically, S-shaped 
as well as asymmetric warps have been shown to arise via 
bending instabilities in the N-body studies of Revaz & 
Pfennigcr (2004); also these have been proposed to arise 
due to gas infall (Lopez-Corredoira et al. 2002). 

In this paper we propose and study the generation of 
the asymmetric warps in disk galaxies. We propose that an 
asymmetric warp arises due to a dynamical wave interfer- 
ence between the first two stable bending modes(i.e. m = 
and TO = 1) in a gravitating system. We do not deal with 
the orign of these modes, we assume that some excitation 
mechanism - perhaps an external perturber, extragalactic 
magnetic field, or gas accretion, or internal bending in- 
stabilities - gives rise to both these modes and hence we 
do not include explicitly any source term in our analysis. 
It is reasonable that, once excited, these modes will in- 
teract with each other. We study a general case so that 
the modes need not even be triggered in a commensurate 
fashion. A different epoch of onset for the two modes for 
the two cases will merely be reflected in a different initial 
ratio of the amplitudes of the modes. 

It was demonstrated numerically by SC88 and Sparke 
(1995) respectively that warping mode {m = 1) and 
the bowl-shaped mode (m = 0) are stable in a cold, 
thin self-gravitatiing axisymmetric disk. The basic idea 
of the present paper is to use this fact and then let 
these modes interfere with each other, and see if this re- 
sults in an asymmetric warp. This idea of superposition 
of modes to explain the origin of asymmetric warps was 
first proposed by Sparke (1995), and also more recently by 
Lopez-Corrcdoira et al (2002) and Castro- Rodriguez et al. 
(2002), but they did not work out further details. 

We do not address here the other important aspect 
of the warping phenomenon, namely its maintenance, 
which is a long-standing problem. In the modal approach 
some reasons for the reduction of warp lifetime are: back- 
reaction of the uniformly processing disk onto a live dark 
matter halo (Nelson & Tremaine 1995), and disk thickness 
and random motion (Sellwood 1996), which can damp out 
the disk warping in a few orbital time scales. Some numeri- 
cal simulation studies (Dubinsky & Kuijken 1995; Binney 
et al. 1998) have confirmed that the dynamical friction 



due to the oblate dark matter halo damps out a symmet- 
ric warp. On the other hand, the maintenance problem 
does not arise if there is gas accretion of ~ few MQyr~^ 
repeatedly generating warps as discussed above. 

In this paper, we have isolated and studied one well- 
defined physical aspect of the problem, namely that if the 
m=l and m=0 modes are excited due to any of the mech- 
anisms mentioned above- namely, tidal interaction, or gas 
accretion, or bending instabilities, then it is inevitable 
that these will get superposed and we study the result- 
ing behavior. Further, even if these modes are short-lived, 
if these can be excited repeatedly, then we can explain 
why a large fraction of spiral galaxies show asymmetric 
warps. 

In section 2, we formulate the equations for to = and 
TO = 1 cases, and give the numerical scheme for solving 
these, and also specify the input parameters for the disk 
and the halo. The results are described in Section 3, and 
Section 4 summarizes the conclusions from this paper. 

2. Formulation: 

2.1. Dynamics of bending modes: 

We consider various possible bending modes in a cold, 
self-gravitating, thin axisymmetric disk in the presence of 
a rigid or non-responsive dark matter halo. The disk,with 
a radial surface density profile E(r), rotates in the equa- 
torial plane {z — 0) of the spheroidal dark matter halo 
with angular speed fl{r) about the halo's symmetry axis 
(r = 0). Here {r,(p,z) are the circular cylindrical polar 
coordinates. The dynamical system described above is the 
same as in SC88 and most of the formulation follows SC88. 

The dynamical equation of a small bending of the disk 
perpendicular to its unperturbed plane(2; = 0) is given by 

^+f^(r-)^J Z = J'self+^halo (1) 

where the small bending is described by a single function 
Z{r,(p,t) and J^seif is the vertical force due to the bent 
disk itself. Thaio is the vertical restoring force near the 
disk plane {z = 0) due to the dark matter halo. 

We consider the small bending of the disk normal to its 
plane as resulting from the linear superposition of differ- 
ent bending modes each described by an azimuthal wave 
number to. So the linear superposition allows us to write 

Z{r,(p,t) = ^A^Zm,{r,y^,t) (2a) 

m 

In the above relation is a constant number which 
denotes the coefScients of superposition and 

Zm{r,v,t) = 3?{/i™(r)e'(--*-™^)} (26) 

describes the bending mode of the disk with azimuthal 

wave number to, where hm{r) is the unmodulated ampli- 
tude of the m*'' mode. Thus, the total amplitude of the 
TO*'' mode in the superposition is given by hmir)- 
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The general dynamical equation for the bending mode 

with the azimuthal wavenumbcr m is then obtained by 
substituting the form given by eq. (2b) into eq. (1) : 

{wm - mfl{r)f - Vh^ir) hm{r) = 



POO 

G / W)H{r, r')[hm{r) - /i„(r')]r'dr' 

f-OC 

+G E{r')Kmir,r')hm{rydr' {3a) 
Jo 



The angular speed fl{r) appearing in above equation gets 

a contribution from the warped disk as well as from the 
dark matter halo described in sec. 2.3. The definitions of 
H{r,r') and Km{r,r') are as follows: 



Km{r,r') 



,t') = / 
Jo 



[r^ + r' — 2rr'cos V' + zq^] 2 



■dip (36) 



/ " '"'^"^^^ ^d^P (3c) 

[r^ + r' — 2rr'cos-0 + zq^] 2 



These functions can be further expressed in terms of com- 
plete elliptical integrals. Note that for the bowl-shaped 
mode(m = 0), the function Ko{r^r') vanishes. 

Eqn(3a) is obtained by substituting in eqn(l) the fol- 
lowing: 

^Zioir,V,t) = -Vh^{r)z„^{r,ip,t) (4a) 

where is the vertical frequency due to the dark matter 
halo in the unperturbed disk plane. 

J^Zif{r,'P,t)=-Gj^ i:{r'ydr' 

'■2" [z„(r,</j,i) -z„(r',(^',t)] 



-V (46) 



[r^ + r'^ — 2rr' cos((/3 — (/?') + Zq^ 

In the above equations zq is the softening parameter in- 
troduced to make the integrals regular at r = r' . zq 
can be interpreted as a finite thickness of the disk. The 
above integral eq.(3a) can be solved by recasting it into 
a matrix-eigenvalue problem. On a uniform grid with N 
radial points the resulting eigenvalue equation for the m*'' 
mode takes the form 

N 

[ujyn^ -2mQ.{ri)u:rn]hm{ri) = ^kijhm{rj) (5a) 
Where 

Kij = Mij + Sij (uh^irj) + Vdc^{rj) - m''^{rjf'^ (56) 
M,j = ArGS(rj) [H{ru Tj) - i^™(r„ r^)] (5c) 



And 



The above eq. (5a) can be rewritten in a more compact 
form: 



{Wrn^I + UmR + A) /im = 



(6) 



Where I, R and A arc the three NxN real square matrices. 
hm is the eigenvector corresponding to the eigenvalue ujm- 
The matrix elements are = Sij, Rij = —2mfl{rj)6ij 
and Ajj as shown above in eq.(5b). 



2.2. Asymmetric warp 

The bowl-shaped mode is an axisymmetric bending of the 
disk and the integral-sign warping mode is anti-symmetric 
in the azimuthal angle (ip). In the linear regime each of the 
two bending modes behaves much like an independent os- 
cillator in the system, so that we ignore any kind of energy 
cascading. This is shown to be valid in the last paragraph 
of sec. 3. By making a linear combination with proper 
time dependence of the two bending modes we are able 
to produce an asymmetric bending mode of the initially 
assumed axisymmetric disk. Thus according to eq.(2a) the 
asymmetric bending mode is described by 



i{r,ip,t) = AoZo{r,ip,t) + Aizi{r,ip,t) 



= Ao^{ho{r)e''^°'} + Ai^{hi{r)e 



i{ujit—ip) 



} 



So that at any time t = r wc can analyse the be- 
haviour of the resulting asymmetry. It was proved that 
the TO = mode and the to = 1 mode in a disk described 
above are stable (sec Hunter & Toomrc 1969), and their 
stability was further confirmed by numerical work (see 
SC88, Sparke 1995). Our numerical work also confirms 
the earlier result. Thus the resultant asymmetric bend- 
ing mode can be assumed to be stable and it establishes 
a classic example of dynamical wave interference pattern 
in a gravitating medium. The degree of asymmetry in the 
disk bending depends on the amplitude of the bowl-shaped 
mode which can be controlled by the free parameter Ao 
in the problem. Thus by varying Ao one can in principle 
produce a rich class of asymmetries that are seen in the 
observations. Recent observations show that asymmetric 
warps arc indeed common and they exist in poor envi- 
ronments too (Garcia- Ruiz etal. 2002). This suggests that 
tidal interactions are not the only process that can pro- 
duce an asymmetric warp. With our approach, we show 
how an asymmetric warp can be generated purely from 
the internal disk dynamics- assuming that these modes 
are generated due to some mechanism. 

2.3. Disk-Halo system 

We have adopted a truncated exponential disk with a cos^ 
tapering as introduced by SC88: 



poo 

VdM) = G i:{r')Hin,r'ydr' {5d) 
Jo 



So e 



-r/Rd 



= r > rout 



r < Ttr 
COS {- 



rtr 



2 Tout Ttr 



-} nr 

I 

(8) 



< r < r„ 
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where is the exponential disk scalelength and S is the 
disk surface density. The gradual truncation of the disk 
avoids spurious modes (SC88). In the calculations (Section 
3), we have taken rout and rtmn to be fixed at 6 and 5 disk 
scalelcngths respectively. 

The flattened dark matter halo is taken as a screened 
isothermal halo which gives an asymptotically flat rotation 
curve (de Zeeuw & Pfenniger 1988) : 

p(r,z) = (9) 

where Rc is the core radius and q is the halo flattening 
parameter. 

The terminal velocity of dark matter halo is given by: 
V, = AnGp.R^V-^^^^^^^ (10) 

2.4. Numerical details 

Equation (3a) is a linear integral equation which is equiv- 
alent to an infinite dimensional eigenvalue equation (6). 
To simplify this problem, a galactic disk of finite ra- 
dius (Vout) is approximated as a system of N uniformly 
spaced rings. Then the resulting problem reduces to an 
N-dimensional eigenvalue problem in the eigenvalue Um 
and N-dimensional eigenvector hm describing the shape 
of the bent disk. This requires the solution of an NxN 
matrix. The softening parameter Zo is introduced to make 
the self-gravitating integral (eq.[4b]) regular along the di- 
agonal line of the N-dimcnsional matrix. We have carried 
out the numerical computation for different values of the 
matrix dimension N. We note that the eigen frequency 
of the ground state of either mode (m = & rr?, = 1) 
is not very sensitive to the matrix dimension N or the 
softening parameter Zo in the problem. Numerical results 
show that as the value of N increases or the softening 
Zo decreases, the ground-state eigen-frequency oji of the 
warping mode (m = 1) reduces extremely slowly. For ex- 
ample over an augmentation of 200 rings the change in 
oJi ~ 5.0%. For convenience we have taken N=80 for our 
calculations. A similar calculation was done for the case 
of the bowl-shaped mode (m = 0) but it did not show any 
definite trend unlike the m = 1 case. But the change in 
the frequency(a;o) occurs only in the third decimal place 
and its quite arbitrary. The insensitivy of the ground-state 
mode frequencies (wo & wi) on the softening is obvious be- 
cause these are the integral properties of the disk. 
However the ground-state mode shape, h„i, in either case 
is quasi-sensitive to the softening parameter. A smooth 
bending of the disk was obtained for typical values of 
the softening Zo ^ inter ring spacing. We note that a 
large value of the softening reduces the value of the self- 
gravitating integral, hence it becomes inefficient to make 
all the rings precess in a synchronized manner, thus any 
resultant warping is dissolved. On the other hand if Zo = 
it is obvious that the self-gravitating integral diverges nu- 



merically. Since we are not interested in investigating the 

structure of the either mode in the resolution ^ inter-ring 
separation, we keep the softening parameter Zo ~ inter- 
ring spacing and this gives us a satisfactory result for the 
mode shape. The sensitivity of the mode shape on the 
softening is understandable because a global mode shape 
depends on the local details of the disk. 

For m = 0, i.e. bowl-shaped mode, the eigenvalue 
equation (6) is linear in Wo^- But for m = 1, i.e. a warping 
mode, the problem reduces to a quadratic eigenvalue prob- 
lem (known as QEP in the literature) in lui . The standard 
way to solve the QEP (6) for the m=l mode is to reduce 
it to a generalized eigenvalue problem (GEP) of the form 
Ax = ojmBx of twice the matrix dimension 2N see Bai et 
al. (2000). This GEP is commonly called a linearization of 
the QEP (6). We solve the linearized QEP numerically by 
using standard technique for diagonalization. For further 
details about solving QEP the readers are referred to Guo, 
C. -H (2004) and Higham & Kim, (2000). 

3. Results 

Recent observations have shown that the galactic disks 
bend with a variety of shapes starting from symmetric 'S' 
shaped warps to asymmetric warps and also 'L' shaped 
and 'U' shaped warps (Section 1). We next show that this 
rich variety of warp shapes observed can be explained nat- 
urally right from the internal disk dynamics by a super- 
position of the m = 1 and m = modes. The superpo- 
sition picture proposed and studied here is verified for a 
few cases in the results from the N-body simulations, and 
there is evidence for an asymmetric warp arising due to a 
superposition of m=l and m=0, or m=l and m=2 modes 
(Revaz & Pfenniger 2005, personal communication), this 
requires further study using N-body simulations. Since we 
superpose two modes, the range of values spanned cover 
the two frequencies wq, wi, the superposition amplitudes 
A{), Al, and the two phases. The unmodulated amplitudes 
hi and ho (see eq. [2b]) are free upto a linear multiplica- 
tion factor. These and the frequencies are obtained from 
our modal analysis. The ratios of ujo/loi are seen to vary 
from 8-20 when the flattening of the dark matter halo (q) 
varies from 0.5 to 0.9. (see the discussion at the end of 
this section) . We do not have any apriori physical basis to 
choose the values of Aq and Ai, these are set by the gen- 
eration process of the modes- hence we cover a reasonable 
range for these. 

All the numerical calculations are done in units of 
G=Md=i?d=l where is the disk mass. In order to show 
the basic results first, we have fixed the dark matter halo 
parameters at Rc = 2.0 ( in units of Rd) , Po = 0.0117 
( in units of M^/R^) with a halo flattening parameter 
q = 0.7. These give rise to the terminal velocity of the 
halo as Vt = 0.68 ( in units of y/GM^jRd). 
In order to check our numerical calculations with the pre- 
vious works by SC88 and Sparke, 1995 we have reproduced 
the ground state mode shape for the modes m=0 and m=l 
for some typical disk and halo parameters as can be seen 
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in Figure 1. The amplitude of these modes are measured 

with respect to the inner disk plane and these wc call as 
the unmodulated amplitudes. In addition to this, Figure 1 
also shows the sensitivity of the amplitudes of these modes 
with respect to disk edge. To study the sensitivity wc have 
varied rtmn keeping the disk outer edge fixed. The solid, 
dashed and dotted lines are for rtrun=5-0, 5.5 and 5.95 in 
both the figure panel. The disk density starts deviating 
from the exponential from r = rtmn and smoothly tapers 
to zero at r = rout- The dotted curves show the behaviour 
of the modes when the disk is almost abruptly truncated 
and in this case too the two modes m=0 and m=l be- 
have in the same way. Also note that the amplitudes of 
m=0 modes are slightly less than that of m=l modes. The 
eigen frequencies corresponding to these ground states are 
not quite sensitive to the behaviour of disk edge which 
can be seen from the SC formula for the modified tilt 
mode(eqn(21) in SC88). This actually confirms what we 
get numerically for the m=l mode. However for the m=0 
modes the cigcn frequencies vary from 0.1315 to 0.1296 as 
ftrun varies from 5.0 to 5.95. 

In Figure 2 we have shown four subplots of an asymmetric 

warp generated with various values of the controlling pa- 
rameter C, = Ao/Ai which we will use as an indicator of the 
relative strength of the bowl-shaped mode, = 1-0 and 
j4i = 1.0 connotes that the bowl-shaped mode and the 'S' 
shaped warping mode both are present with their full un- 
modulated amplitudes, hm (see eq. [2b]). These plots are 
made for an epoch of To/8 as an illustrative case, where 
To = 2i^/uJo is the characteristic time scale in the problem 
and uJo is the ground state eigenfrequency of the bowl- 
shaped mode. As wc move from Fig. 2a to Fig. 2d, the 
value of the controlling parameter ^ increases, that is, the 
relative strength of the bowl-shaped mode increases. Since 
in our analysis bowl-shaped mode is the only cause of 
asymmetry in the bending, we can see from Fig. 2d, that 
the warp is most asymmetric i.e. the degree of asymmetry 
is the highest. 

In Figure 3, we have produced asymmetric warps at differ- 
ent epochs {t) keeping the controlling parameter Q = 1.0. 
The values of the sampling times t were so chosen that as 
we move from Fig. 3a to Fig. 3d, we see the diverse phe- 
nomena of disk warping including symmetric as well as 
asymmetric ones. This is due to the characteristic oscil- 
lation of the disk in a bowl-shaped mode. The disk oscil- 
lates from 'cupped upward' to flat, to 'cupped downward' 
and back again with its typical ground state frequency uJo ■ 
During this flexing of the disk it gets a chance to interfere 
with the S-shaped mode of the disk and this produces a 
rich class of dynamical asymmetric figures. In a real galaxy 
the bowl-shaped mode may decay due dissipative effects. 
Even then the mild asymmetries in disk warping which 
are seen in observations can be considered as a signature 
of the bowl-shaped mode. 

In the above two figures (Figs. 2-3) we have shown the 

warping of the disk with respect to the inner unwarped 
disk plane because in actual observations the warping is 
quantified based on measurements with respect to the in- 



ner disk plane. Next, we define a quantitative measure of 

the asymmetry seen in a disk galaxy, the asymmetry in- 
dex, as follows (see also Sanchez-Saavedra et al. (2003)): 



a. 



asym 



Wright - O^leftl 
aright + aieft 



if a 



right 



+ a] 



eft 



7^0 



where aright is the angle between the line joing the cen- 
tre to the outer most point on the right hand side of the 
particular warp and the inner disk plane, aieft is defined 
similarly. Thus aaaym is the normalized or fractional value 
of the asymmetry. In terms of this definition, the results in 



Figures 2 and 3 give an asymmetry of aas 



0.2 - 0.4. 



Thus our results for asymmetry agree with the typical 
range of measured values in Sanchez-Saavedra et al. (2003) 
and the asymmetries when measiircd from the figures 
given in appendices (A & B) in Schwarzkopf & Dettmar 
(2001). For example the value of asymmetry index for 
the S-shapcd asymmetric warp(galaxy name- AM 1134- 
323 (R)) when measured from Fig(6) in Schwarzkopf & 
Dettmar (2001) gives a value 0.33. This agrees closely 
with what wc get from Fig(2). Of course true compar- 
isons with observations is a bit tricky because of the free 
parameters(j4o and Ai) involved in our problem. Our ba- 
sic idea is to show that with this modal superposition ap- 
proach, one can reproduce various asymmetric warps that 
are seen in observations. 

Clearly when aright = ai^ft the asymmetric index 
Oiasym vanishes. This implies a purely symmetric warp- 
ing of the disk. Note that eq.(ll) avoids the possibility of 
oiright and aje/t being simultancosly zero as the case is not 
physically meaningful. When aright = or Q;;e/t = we 
obtain unity value of the asymmetric index. So aasym = 1 
(maximum asymmetry) denotes a perfectly 'L' shaped 
warping of the disk. This can be called a one-sided warp. 
The values of aasym vary from 0.0 to 1.0. It is obvious 
that most of the geometrical shapes of galactic warps seen 
in observations can be classified using a single parameter 
aasym- Jimenez- Vicente et al. (1997) used a different set 
of 3 parameters to describe the shape of warps. However 
neither their model nor aasym alone can describe the pe- 
culiar shape of the warp of our Galaxy. But aasym being a 
single parameter is very useful in depicting the underlying 
asymmetries in warps. Given the extreme importance of 
this parameter, we next study the dependence of aasym 
on the dark matter halo parameters. 
Figure 4 shows the variation of the asymmetry index, 
aasym, due to the variation in q, the flattening param- 
eter of the dark matter halo for a fixed value of the con- 
trolling parameter (C=l) and at an epoch t = To/8 as a 
special case. As the flattening parameter of a halo of flxed 
mass is varied its central density, po. and the core radius, 
Rc are bound to change. Therefore we need to calculate 
the Po and Rc as functions of q. By imposing two con- 
straints: the mass within a thin spheroidal shcll(Binney & 
Tremaine 1987, pg 54) and the terminal velocity {vt, see 
eq. [10]) of the halo should be independent of q, we obtain 
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Po(<?) and Rc{q) in terms of their spherical counterparts 
as (see Narayan et al. 2005): 

po{q) = Po{i - i—r) ' Md) = Rc{i)( 

(12) 

where e = (1 — q^Y^^. As the halo becomes more flat- 
tened (smaller q), the resulting asymmetry aasym goes 
down which means that the warp becomes more sym- 
metric. Conversely, as the halo becomes more spherical 
the asymmetric index aasym goes up showing that asym- 
metric warps are more likely to be found in a less oblate 
dark matter halo. This can be explained as follows. As a 
constant-mass halo is flattened, the ground state eigenfre- 
quency (wo) of the bowl-shaped mode increases, so that 
the disk flexes up and down more rapidly. Therefore the 
characteristic timescale (jo) of the mode decreases and 
thereby its unmodulated amplitude reduces faster at the 
measured epoch, this in turn reduces the asymmetry of the 
warps. While this figure illustrates an important physical 
point from the model, its direct observational verification 
is not possible. 

In Figure 5 we plot the variation of the asymmetry in- 
dex, aasymi with the core radius Rc of the dark matter 
halo (given in units of Rd), again for a fixed value of the 
controlling parameter (C=l) and at an epoch t = To/8 
and for a halo flattening oi q = 0.75. In contrast to Fig. 
4, here as Rc increases, the halo mass within a given ra- 
dius goes up. Fig. 5 shows that asymmetry in warps are 
likely to more in a halo with a smaller core radius; while 
as Rc increases, aasym goes down. This is because as Rc 
and hence the halo mass increases, the frequency ujo of the 
bowl-shaped mode is raised, and this results in a smaller 
intrinsic or unmodulated amplitude of the bowl-shaped 
mode. So in absence of an appreciable amplitude of the 
bowl-shaped mode, the resulting asymmetry in the disk 
bending also reduces. This result of decreasing asymme- 
try in warps for higher mass galaxies is in agreement with 
observations (Castro- Rodriguez et al. 2002), which further 
supports our model. 

The variation of the asymmetry index aasym on the halo 
mass can be seen from the Fig. 5. It is explained above 
that as Rc/Rd increases and hence the halo mass increases 
the asymmetry index goes down. While as the halo mass 
decreases the amplitude of warping mode (m = 1) starts 
decreasing and in the limiting case when the dark matter 
halo is absent m = 1 mode turns out to be a trivial tilting 
of the disk. And it is hard to generate m = mode in 
the absence of a dark matter halo. Thus in the limiting 
case when halo mass is tending to zero, the asymmetry 
index starts falling down to zero and however as we have 
checked that this happens when Rc/Rd starts becoming 
a fraction which is not a physically acceptable regime of 
Rc/Rd for any reasonable disk-halo system known. 

So far we have considered a screened isothermal {p oc 
at large radii) dark matter halo, producing a fiat ro- 
tation curve, for the generation of asymmetric warps. The 



flaring of HI in the outer region of our Galaxy favours 
a dark matter halo with a steeply falling density profile 
(p oc or r~'^ at large radii) ( Narayan et al. 2005). 
Note that p oc r~^, at large radii, is the popular NFW 
density profile (Navarro et al. 1996) for the dark matter 
halo. We have checked that the amplitudes of the two 
modes (m = & m = 1) and their ground state eigen 
frequencies behave in a similar fashion in both cases when 
p oc and p oc or at large radii. Hence we can 
say that qualitative behaviour of the asymmetry would 
remain similar. However we would like to carry a detailed 
investigation of the asymmetric warp with respect to dif- 
fernet halo profiles in a future work. 

In Figure 6 we show one of the extreme form of asym- 
metry of warp shapes, namely an L-shaped mode (see 
above discussion following eq.[ll]). In principle this can 
be achieved by reinforcing the net amplitude of bending 
resulting from the superposition of the two modes to zero 
on one side of the disk. We get an approximately L-shaped 
mode for = 0.55 and Ai = 0.20 for the usual input val- 
ues of ^=0.76 and Rc=2.0. Clearly the root of such an ex- 
treme form of asymmetric warp lies in the fine balancing 
between the two modes in the process of superposition. 
Nearly 5 % of the asymmetric warps observed have this 
structure (Sanchez-Saavedra et al. 2003), which can be 
explained by our approach. 

Another form of deviation from an S-shaped distribu- 
tion is a U-shaped warp, which comes about naturally 
when the amplitude Aq of the bowl-shaped (m = 0) 
mode is enhanced w.r.t. the amplitude Ai of the warping 
(m = 1) mode (see Figure 7). Nearly 7 % of the asymmet- 
ric warps observed have this structure (Sanchez-Saavedra 
et al. 2003). A little juggling of Ai and Aq values is re- 
quired to obtain the L-shaped and the U-shaped modes, 
unlike the case of the other asymmetric warps. This per- 
haps explains their lower observed frequency. 

Thus our approach can naturally produce the variety 
of asymmetry seen in disk galaxies, including the pecu- 
liar L-shaped and U-shaped warps. The resulting values 
of asymmetry (Figs. 4-5) agree with the range of observed 
values (Sanchez-Saavedra et al. 2003 and Schwarzkopf & 
Dettmar 2001). 

In order to compare the resulting asymmetry with the 
observed values, it is useful to obtain the true value of tq. 
We have checked that this lies in the range of ~ 0.5 — 1.0 
Gyr and n lies in the range of ~ 6— 12 Gyr, for the realistic 
galaxy models with a range of i?^ 2 — 3 kpc, Mn, the 
disk mass~ 2 - 6 x 10^°Mo, R^ ~ 2Rd, q = 0.7 and the 
ratio of the halo mass to the disk mass is~ 1.4 within the 
outer radius 6Rd)- The ratio of the frequencies luq/loi is 
~ 10, note that this value is obtained for q=0.7. For other 
values of q the ratio ujq/u)\ varies. Numerical calculations 
yield a range of values varying from 8-20 for q varying from 
0.5 to 0.9 respectively. In the absence of any dissipation, 
the superposition of the two modes will show a cyclical 
variation and hence the above values represent the typical 
asymmetries seen. Even if there is dissipation, so long as 
there is a regeneration of the modes (see Section 1), then 
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such asymmetric warps can recur and hence can be seen 
effectively over much larger timescales. We note that since 
the natural frequencies of the two modes are so different, 
we are justified in treating the problem as a simple linear 
superposition rather than having to take account of the 
mode coupling. 

4. Conclusions 

We propose and study the origin of the asymmetric warps 
in spiral galaxies as arising due to the superposition of 
the standard S-shaped warps (to — 1 mode) and a bowl- 
shaped warp (to = mode). This is a simple but general 
model for generating asymmetric warps that is indepen- 
dent of whether these modes arise due to tidal interaction, 
or gas accretion, or bending instabilities, or any other per- 
turber. We do a modal analysis of a disk embedded in a 
dark matter halo, and obtain the solutions for these two 
modes which are then linearly superposed. In the charac- 
teristic oscillation of the to = mode, the disk oscillates 
from 'cupped upwards' to flat, to 'cupped downwards' and 
back again. During this flexing of the disk, it gets a chance 
to interefere with the S-shaped (to = 1) mode of the disk, 
and this produces a rich class of dynamical, asymmetric 
warps. The results obtained naturally explain the wide 
variety of asymmetric warps observed in spiral galaxies, 
including the peculiar L-shaped and U-shaped warps. 
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Figure 1. The ground state mode shapes and their sensitivity 
w.r.t. disk edge for the two modes rn=0 and m—1. The mode 
amphtudes are in arbitrary units. The solid, deished and dotted 
Unes are for rtrun=5.0, 5.5 and 5.95 respectively and rout=6-0- 
Here, the dark matter halo flattening q = 0.7 and the core 
radius Rc = 2.0. 
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Figure 3. A plot of the bending amplitude of the asymmetric 

warp versus radius, at difi'erent epochs t. As t goes through 
the cycle of the bowl-shaped mode it generates various kinds 
of disk warping from an asymmetric to a symmetric and then 
an asymmetric case. Again, the dark matter halo flattening 
q = 0.7 and the core radius Rc = 2.0. 



a 



asym 




0.5 0.6 0.7 0.8 0. 



Halo flattening q 



Figure 2. A plot of the bending amplitude of asymmetric 
warp versus radius, for various values of the parameter ^ = 

Ao/Ai denoting the ratio of the superposition amplitudes, at 
an epoch t = To/8. As ( increases to 1, the asymmetry becomes 
prominent. licrc, the dark matter halo flattening q = 0.7 and 
the core radius Rc = 2.0. 



Figure 4. A plot of the asymmetry index, aasym, versus q, 
the dark matter halo flattening. As the halo becomes more 

spherical, the asymmetry in warps goes up. Here, the dark 
matter halo is flattened in such a way as to conserve the ratio 
of halo mass to disk mass within optical radius of the disk, and 
the terminal velocity of the halo. 
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Figure 5. A plot of the asymmetry index, aa 



versus Rc 



the dark matter core radius. This sliows that as the halo core 
radius and the halo mass increase, the asymmetry in warps 
reduces. Here, the dark matter halo flattening is kept constant 
at g = 0.75. 



Figure 7. A plot of bending amplitude versus radius which 
shows an U-shapcd asymmetric warp produced by sotting /lo = 
1.0 and^i = 0.1 at an epoch t = To/8. The dark matter halo 
flattening is kept at g = 0.76 and Rc = 2.0. 




Figure 6. A plot of bending amplitude versus radius which 

shows an L-shapcd asymmetric warp produced, by setting 
^0 = 0.55 and Ai = 0.20 at an epoch t = To/8. The dark 
matter halo flattening is kept at g = 0.76 and Rc = 2.0. 



